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SYt^OPSIS 


This thesis <ieels with a general mxn Job-Shop 
3chec;ulirig probleE vvitli an objective to minimise the 
sum of the vcosts of lob late— times' and machine - 
i die- times. In optimei procedure' u sing the branch 
Ohd bound technique in which the tree of sequences 
is constructed on the basis of clashes in the 
Gnatt-chert of tire problem has been evolved. Also 
S'.ibs>ptiraal procedure., based on tire proposed 
tocirniouo, has been suggested to replace tire thumb 
and irfc'ur?istic rules used hitherto. 
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CiiiiPTEfi I 


THE JOB-SfiOP SCHEDOLI^jG PROBLEM 

In recent ye.- rs tt'-ere has been s great interest 
in Industri'c.1 scheduling problems from the point of viev^ 
ol both practice ann research. Scheduling is the 
managerial decision of fixing the sequence of some 
activity (job) on some facility (machine) v.ith respect 
to otiuir activities (jobs) optimising certain cost or 
profit function. In this introauctory chapter an attempt 
i s made to enumerate the vast dimensions of the above 
p roblem. 

A job- shop in any industry usually consists of 
folloVving ; 

(1) A number of various kinds of machines, each kind 
of machine may be avr.ilcible in more thutn one 
number. 

(2) Ihe rrvi material, i tsastorage, arrival and 
distribution mechanisms. 

(?) Hie human resources and their efficiencies, 
w) ■ The storage and despatching mechanisms for the 
finished products. 

The jobs, to be processed in the above job- shop 
usually have the following sp scifications : 

(1) The sequence of machines to be used by the jobs. 



(2) Tb.e ^orocesjjing times of jobs on various machines. 

(3) Cost function of the completion of the jobs. 

(4) i'rrivel pettern of jobs (or job orders). 

The object of the floor manager is to manage the 
job- shop most econoraically. Inclusion of all the above 
factors and r. fevt more ma'kes a problem of managing the 
complete industry. Houever, tne job- shop-problem, as it 
is generally referred to, is that of assigning jobs on 
various machines. Before describing the job-shop problem 
in the present context in gre.- ter detail, let us see the 
kind of the cost functions associated vvitlr this problem. 

The cost functions are the result of the inter- 
action of the shop nuno.gement and the customers. Following 
may be tho usura customer’s das ends : 

(1) The custcr..er asks for the finished job on a due date. 
If the job is late the penalty as e function of the 
delay is specified. This function may be nonlinear. 
Efrly dispa^tch may or nay not be favoured. 

(2) The customer wemts the job at a certain date, what- 
ever tho cost be. If it cannot be given he might, 
like to scrap the deal. No ea^rly despatch may also 
be specified. 

( 3 ) customer wants his job to cost least no matter 
when it is completed. 

(4) Another customer wants his job to be rushed. They 
may be prep.' red to reward early aespatch. 



Mf^ny other sitiu^tions or combinetions of them 
cm occur. Being ff ced vvith varying demsnds like above^, 
the shop mccnrgement has to sssess their capacity. ihey 
hrve to prepare the tentative schedules Vvhich together 
■with the above demands also incorporates the cost of 
machine idle times, and inprocess job inventories. After 
accounting for the allovvances for mishaps, failures etc. 
the management should be a.ble to say ^h.ether it can make 
them or not. If he cen, then he should be able to 
specify the costs incurred and have a diologue with the 
customer. i\ll the above cost requirements are nonlinear 
and therefore requires nonlinear programming techniques 
for getting an optimal solution. 

Besides tire complex nonlinear optimizing function 
the constraints on the job— shop a.lso pose a serious 
problem. Let us enumerate them as follovjs : 

(1) jNo machine may process more then one operation at 
8 time. 

(2} Certain operations on the jobs may be interruptible 
and others not. 

(3) Assembly or diasscmbly of some jobs. 

(4) Dynamic or batch input of jobs. 

(5) Time bounds for the gap between tm operations 
of a job like the usual examples in chemical 
processing. 

(6) Oncertainty in the processing times given by the 
customer, or estimated elsewhere. 



(7) One or raoro machines of one kind. 

(8) iJ-tcrnrte sequence of processing on jobs, i.e., 
oltornete routing of jobs. 

(9) /. job nry be abendoned if oraer is cancelled. 

(10) i', job Dicy or ra&y not be processed by more then one 
machine. 

(11) Some jobs may not hove rigid routings. 

The list is not exhaustice but only illustrative 
of hovv complex: the practical problems can be. There 
exists no analytical model to consider the above 
mentioned general job- shop scheduling. Even the 
simplest problems consisting of five to six machines and 
the seme order of jobs require n large computer time to 
obtain idle optimal solution. only recourse left is 

to use heuristic rules and simul^i-® the job-shop. It 
requires considerable experimentation. 

In this ti^esis an optimal procedure, based on 
the branch and bound technique, has been developed for 
the solution of a general job-shop scheduling. The 
job-shop considered is as follo'^'''® : 

(1) Batch input of jobs. 

(2) i\ll machines are of different kinds. 

(3) The processing times are known apriori. 

(4) The costs are linear functions of jobl&teness 
and txhe machine- idle times. 


(5) el tern 'te routes. 

(6) I'Jo riesc.-nbly OT d-tasseably and no interruptions of 
er.y op orations. 

L^ter a sub-optimtl procedure is given based on 
true rbavo procedure. It c‘-n be substituted for any 
heuristic rules of scheduling in use so far. In the 
l-'ust sorie suggestions r re given for taking care of the 
riOiillnG”r cost functions, i.e., due dates, interruption 
and './ore tron crie nocbine of one kind. 



Cfi/’.PTER II 


PRESENT A4ETH0DS OF SOLO TIDE 

Kiio'...ing the complex structure of a general job 

shop probleiTi;, no one has tried to solve it completely. 

Many simplified models have been put forvvard and analysed. 

Freouently considered simplification are listed belovv : 

(1) No machine may process more then one operation at 
a tirae. 

(2) No operation on any job is interruptible and no job 
requires processing by a machine more than once. 

(3) No o.ssoTnbiy or c.i assembly csllovved. 

(4) The processing times include the set up times and 
'•re fixed. 

(5) There are no constraints on the gap between two 
operations on a job. 

(6) No alternate routes. 

(7) There is only one of each type of machine. 

(8) There is batch input of jobs. 

The optimizing function has also been simplified, 

and tile ones frequently considered are : 

(a) Minimizing the make span time or the time of 
completion of all processes on trie jobs. 

(b) Minimizing the idle time of the machines. 

The less frequently considered optimizing 

functions are : 

(c) The cost of lateness of jobs and idle time on 
machines together. 



(.c;} I'li'eeting tn'--- ciue ac. t'js oi tr^e jobs. 

Besides considering the simpler problerris, the 
size of the problems considered rre slso small. For s 
m jobs n machines tnere exist (mJ number of possible 
combinf tions. Ihis number is a.round a million for a 
six by six problem. Both anelyticcd and heuristic hti 

rules have been deployed to find optima-l and sub-optimal,- 

aciutior s, vi'ich are revievveu. in folio vuing argicles. 

2 . 1 BOKTS TECfihlQOES 

Despite tire vestness of the problem Giffler et al(l) 
report on coerpu tf tional experience with the complete enumer- 
' tion of all ie.asible scnedules. According to them it 
need se? rcli onl;/ the a’ctive feasible schedules which £;re 
defined as follows : 

In active schedule is one in which bo tie job and 
n': a ; c i line a r a n o t s i mu 1 1 ■ n eo u si y i dl e . 

t'.jch scnedules c-.n be generated by assigning a 
job to a .'uachane at an earliest time at which both the 
j o ]■) n u th ; m a chi n e c- r o f ree. 

It reo.uccs considerably tiie number of feasible 
schoauies to be vsearciaed. hote that the active scheduling 
is dependent- on the distribution of the processing times 
and routings on different jobs. Hence it cannot be 
calculated in advance. Giffler suggests Monte Carlo 
sampling of. the Active feasible schedules. In this 



tochnio’je queues ere formed on the rnechines. If conflict 
occurs ct eny machine it is resolvea by rrndom choice. 

They observed thr t the short schedules hf::ve higher probr- 
bllity in the Monte Ceorlo distribution then in the true 
distribution obtained by linear complete enumeration 
rlgorithms. Ti-e reason propounded by Giffler et si. is 
thrt short schedules have, on the average fevuer conflicts 
to be resolved then long ones. Since in the Monte Carlo 
rlgorlti'im each conflict is resolved et random and therefore, 
the probrbility of obtaining a particular scneoule is 
approximately proportional to the number of resolved 
conflicts it has. Hence the longer schedules are less 
pro!) bio. ./'.fter experimentation with a. large number of 
problems they dreVv time following conclusions i 

1) A Monte Cc'rlo process th....t uses soirie heuristic 
rules as guides in its random choice, v.ill be considerably 
superior to a purely random choice. The “snortest imminent 
operation'' (SIO) has been mentioned to be a, good choice 
for such a rule. In tnis rule tlie operator assigns that 
operation to tne smachine v.hich will be completed first. 

2) Hie larger the ratio of the longest to shortest 
30 b time (tile shortest job time being 1 ), the largar the 
combinatorial aspects oftne problem and therefore more 
difficult it is to schedule. 
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^■'•2 SI:-!'PLE LOADII^G i-MD hEURI^TlC ROLES 

i"rbitr?ry loading rules have so far been the only 
ffiothods in usf- in the industrial scheduling ’situations. 
FolloV'.'ing rrc: some of these rules to pick up the next job 
fron £ ououo for sch. ecu ling : 

(1) First come first sor\''e 

(2) Job with the shortest completion time 

(3) Job tir^ct hc'.s the maximuia number of uncompleted 
operr.itions still to be performed, etc. 

Is against these rules, some rules are devised 
based on tlie skill and the experience of tiic schedulers. 

Tlrn so are termed heuristic rules. Some of them ere as 
folio -s : 

(1) SIO :: Choose that oper-ation for next scheduling 

^ivhich cm be ccmpleted in shortest time. 

(2) LET : Choose tirc^.t job for next scheduling which 

has longest remaining time for completion. 

(3) JOB SLACK : Ctioose that job for the next scheduling 

which had so far ttie largest ratio of slf?ck time 
to operation. This ratio is calcul cited by 
dividing tire slack time on a job by tiie number 
of operations so f'-r scheduled on that job. 

(4) Expected delay time : In this due-date problem the 

total delay in any job is estimatea heuris- 
tically on the basis of the delay already 
occured on that job. The job vilth the largest 



delsy time end its essocirteci cost;, is chosen 
to be scheaulea next. 

(5) JOB SLi^CR Ratio : It is calculated by aividing job 

slrf.cb hours by iiours remaining until the due date. The 
V. Itj'lrrg^rSt ra.tio is chosen for next scheduling. 

V. Gere (2) i-ps studied the performances of ell the 
rbovc rules together Vvith random selections rule for the 
nriy e-date'' problem. After selecting an operation according 
to some rule, he uses th'e follovving enhancement schemes : 

1) j.itorni' te Operation : ... Schedule -the operation 
according to the rule. Kov. check to see if this makes 
another job 'ctiTical' (tti: t is to sec if the slcock of any 
other job has just becomo negative, or if positive, has 
re.-'cnrd ccitcin critical leval). If so, revoke the last 
opjcration ''nd scl’..'a''cjule thao next operation on tice cbitical 'job. 
Check ag-in for la.tcness. If scheduling the second operation 
does not cause any job to be critical, whereas the first one 
did, tJuen schedule tire second one, otherwise schedule the 
firarb one(.tl..at is one vvhich was aictated by ti"-e rule). 

2) Look ohecd : ... vlaen an operation is scheduled, 
v.c may ask : Is tiiere a critical (i.e. late or nearly late) 
job due to reach this machine at some near future i. our, yet 
before the scheduled oporc-tion is completed? If so, schedule 
that job. ctieck to see the effect of this on other jobs. 
Eitbier let the schedule remain » or replace the opera-.tion 
v/ith tile previously scheduled opera tion, depending upon 

the resulting job lateness effected ... 



3) Insert : Once s ’’look ahead” Job has b-een 
scheduled there is a period of idle time on the machine, 
starting at the 'present time’. If there is a Job in 
the queue whose next operation can be completed by the 
time the look ahead job is due to &.rrive at this machine, 
then the open' tiori. should be scheduled. . . . 

Gere's results arc vt-ry interesting : 

(1) The selection of e priority rule for discriminating 
betVveen jobs competing for time on' the machines, is 
not as important as the selection of a set of proper 
heuristics (i,e., the enlia<n cement schemes). 

(2) ' .The SIO and job slack ratio are the simplest end 

most effective rules. 

The heuristic enhancement schemes proposed by Gere 
-re complex. Ho’wever it is unavoidable because of the 
combinatorial nature of the problem and the nonlinearity 
of tile cost function. 

There is a large list of authors Vvho proposed 
different heuristic rules. Efficiency 'Of these rules depend 
upon particular kind of the problems they have vvorked on. 

The heuristic rules do give quick solutions. 

Eowc-ver, for different optimising functions, different 
heuristic rules have to be invented. No author has 
discussed some generalised procedure to develop heuristic 
rules for different optimising functions. 
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hi c.ttenipt in the fbove direction vkas made by the 
cuthor of ttiis thosiG. The scheduling of an operation can 
be decided by en index which is calculated on the basis 
of the pcremet:; rs ; 

(1) overage tix^e of oper tion on a job. 


(f ) T'^urobe i' of ronic inirg operations on a job. 

(3) The drir.'tion of the ^present’ operation being considered 
for the uecision of scheduling. 


f so.- rcn ur.:s maae by giving aiiferent v^eights to 
these paramet; r5j- The^ weight thus obtained for different 
costs on th.- L;:me prohlou.: anci uifferent problems vvere 
compared. Ih-ei-e soer-ed to bo least correlation tu-nong the 
r-v;sults. Th.e reason for this is surely tiie combinatorial 
structure of ti:.e problem. Tnis attempt vvas abondoned in 
favour of tee hope to be able to find such a heuristic 
procedure from the Branch and Bound technique vhich assures 
an optimal solution. This portion is dealt in the fourth- 
chapter. 


2.5 LEARNING TECHMIQUES 

The cutl'iors working on the heuristic procedures 
of solution ' found that some rule works better in one case 
tham others, in another case tlie reverse may happen. This 
prompted Fisher and Thorxjpson (3) to use more than one such- 
rule to construct a. schedule. 



They have usedilonte Carlo random scheduling by 
varying the probc'bili ties of selecting specific decision 
rules at any stage of the scheduling process. The vari- 
ation is depo-ndent upon the ill story of cctions and results 
encountered by previous decisions at the sarae point. The 
actions are the rule choices^ the result is the maJce span 
time, the time of completion of the v.,hole set of jobs. The 
probability of selecting C'. rule is varied on the basis of 
the comparison betviecn the make span time of the present 
schedule v-dtt' r--spect to tire previous schedule or some 
previously established best scheaule. Fisher et.al. have 
used blO and LftT rules for tneir experimen tstlon. They 
used many different heuristic oevices for improving upon 
tiie probabilities of selection of any rule. They conclude 
that learning is possible and tire learnt combinations of 
tv.o or more rules gives a schedule better than the best tirat 
can be generated using the rules singly. They iieve also 
observed that using SIO more frequently in the beginning 
and LRT in the end of the scheduling process is b . good 
combination. 

It ha.s the effect of getting machines to work 
quickly in the beginning and later concentraiting on the 
long jobs. 

Both the amount and the speed of learning depends 
strongly upon our ability to find good heuristic rules to 
start with. So here agadn the problem comes to finding 
rules suited to the optimising function of tiae problem. 



An altorric'tiv-e to tb.e learning using different 
rules es choices, c 'n bo -tbo use of a vveighted index as 
raentioned in the previous article. The learning can be 
done by assigning probobinistic weights to various pera- 
ceters. The boltleneel with -all types of learning is 
that it bias to he done tfresii on each problem. 

2 , 4 BfiAhai flip BOUrjp TSCfiJjlOUES 

'H'is is a SG'-Toh technioue vnich aims at finding 
a T-^ay (or ways) for the early discarding of some classes 
of prrmu tatlons which arc suspect of being ’bad’. It 
has ]?eon used to find the iiinirrum malre spe'n time of a set 
of jobs. 1^50 St of tho vvort In this direction is limited to 
thio jproblems in Vvhich all the jobs have identical sequence 
of processing on various n cc-rlnes. 

Lomnicli and Broun (4) and (5) have used this 
technique to find the rrdnirnara make span time for the type 
of problem nientioned above- Hey have used hoy’s result 
'which says that ''in searciiing for optimum solution (to the 
make span problem) vee cam limit ourselves to the investi- 
gation of those schedules vlid have an identical order of 
jobs on the first tVKO machines and an identical order of 
jobs on the last t'wo machines. ” It reduces the number of 
permutations from (nJ)’'"^ to (nJ)°^ ^ where n and m are 
respectively the number of jobs and machines. For the, 
3-machine case LomnicM has given the folio wing procedure; 
(1) Calculate the bounds on the makre span by considering 
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each of the given ;jobs to be scheduled first. These 
bounds represents th© nodes. Choose the job- with the 
lovt’est bounc and mark it. 

(2) Calculate bounds for deciding the next position from 

tile remaining jobs. The lowest bound of these and tiie 
unmarked nodes of all previous steps, is chosen anu 
marked. Further decisions are taken from this node 
onward. 


This procedure is continued until lowest of all the 
unmarked nodes is found, which gives the complete schedule. 

Ability of getting a good lower bounds is the crux 
of tine success of this method. Brown and Lomnicki propose 
the following bounds : 

"Gj : total time of processing of all the jobs on 
machine j plus the shortest time of processing a job 
(say ith job) on machines 1,2 . . . (j^l) plus the shortest 
time of processing a job (say 1th job) on machines 
(j 1), ... m, where clearly the job i must be different 
from job 1 


For the tin re e-machines case Lomnicki have used yet 



another result from Roy’s thesis that reversing the order 
of machines has as a solution; a permutation in which the 


of inverted problem 


between the 


duration the more 
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a:L t i'l cult it is, on U:-; ■ v sr-' go, tc evcluote the optim.:.'! 
seouencing. Tt-is is tl-.e probier of bccktracking in the 
solution tree. 

In their peper they hove not si'own ho’vv the generc- 
lisrtion to r raschiious v»ill give the optimal result. From 
Roy’s first conclusion the number of permutotions to be 

... __ o 

secr-ci:.'.(v: ..:r-'' (-U.i)'“ % vi/herens tt:e Obti'iors have investi- 
grted only nj sohecules. The remaining permutations are 
these in ohicn the orders of processing of jobs on 
int ■rrr'Cj.rtc mocL.ines ri'-; .not identical. Later Ignall and 
Lehr'' ge '.'.'.'S vi'cn r.icnod an oxempie to this effect, hovvever, 
he fcolls tiict restricting ourselves to identical orders 
v>ill give. Cgoov., schedule". Si.ich permutations cannot be 
rccount-''c if we take a job as an entity for the scheduling 
decision, he must find a branch and bound procedure which 
is based on scheduling the individual operadhons. 

Ignall and Scharge(6> haive given a branch and bound 
procedure very similar to tiie one given by .Brown ajid 
Lomnicki except thr t certain branches are closed for further 
downwa..rd branching by periorming dominance checks. These 
checks aire used 'when ba cktra-cking becomes necessary. Ttiey 
argue that if two partial sequences involve iaenticsl jobs, 
then the seouence vvith tire larger bound need not be 
considered anymore. Later uudek end 3nitn, Bagga and 
Chakrabarti (17) and G.B.Mcmohen (16) have improved upon 
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theyo '.. orrl nrijce checl^s to be opplieci to tvio ptrtiol 
S(N-'u (riic'/r- ill v-hicli only e. fovv scheaulcd jobs £re in 
conoon. Ironicelly tocsc cliecks depend very strongly 
on the idon tleol sequoricGS of processing on sll jobs. 

Iiiotlv end Tei'J.ton (7) erne out with enoitiher forau- 
Irticn t'u t Korn si' .hos riven a counter example for it. 

B. S. polmer (f ) lies given a near optimal procedure of 
getting on optimor. sequence by c slope-order index of 
(li scri' dnr tion . Tl-e index, if converted in non-math e- 
Tiirtlccl Icngufge, says that "give priority to the items 
having tl e strongest tendency to progress from short 
times to long times in tlie sequence of process". Pedmer 
has liniitoc his technique to problems of finding make 
span timelri vvMch all the; jobs nave identicrl orders of 
processing. 

Some improvements in time lower bounds as given by 
Brown and Lomnicki and otriers, have been proposed by 
J.N.D. Gupta (15) in ills recent paper. He reports 5 to 
30^ gain in reducing the computation time. 


2. 5 INTEGER LIHEAR PROGRMlBG HOLOTlOf^S 

Storey and vtdgnrr has uone the pilot stuay in 
this field. They introduce binary valued variables Xj_^ 
vdiich are given a value 1 if item i is scheduled in 
order— position j and zero otiierwisG. Humber of such 
variaxbles is quite large even for small size problems. 
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They start '.v 1 th sn optlnial simplex solution and add tde 
in tcgan* .cons trains one by one till integer solution 
is found. iliey considered various schemes for selecting 
the basic veriables ; nd its corresponding equation on 
which to .form the aeu cens train t 

(1) Select, c basic variable having the largest fractional 
part. 

(f) Select an closest to 0.5. 

(3) vSelcct an .x. . < 1 closest to 1 . 

■ ±2 

(4) Select a,rj Xj_ (-closest to 0. 

(5} Select an 2:4 ^ 1 closest to 1 ; Vk^henever the objective 

Junction hf-.s a fractional port, force the value of tice 
objective function to tiie next integer level and so on. 

licey observe that rule 3 v^orks bettor than rule 2 . 
Hiey contended tlia.t they have not yet found an integer 
programming method tiaat can be relied upon to solve the 
p ro b 1 em s r a.p i dl y . 

The formulation of a generalised mxn job- shop, 
given by Greenberg (IC), is one of the leading steps in 
obtaining analytical solutions to ticis problem. He has 
used the branch and bouna technique as given by Land and 
Doig (11) to solve the linear Integer programming version 
of the problem.. Going down the tree he fixes the sequence 
between two jobs on a maci’iine. It is equivalent to assign- 
ing a binary vadued variable 2 .ero or one. The bounds are 
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ti.'e GVfl'j otiori of the cost function. /. detc.iled r-nslysis 
of tois procedure is given in Chapter III. 

F.gon ,Bula.s has proposed rn implicit enumeration 
'■■l;,ori thm for solving integer linear programming problem 
(LPP) in vlnich certain branches are bloclied on the basis 
of cert- in tests. The doiratninance checlfs used by Ignall 
and Schrago f. nd their follovn-irs ;'lso do essentie-lly 
the- some thing. The integer LPP formulations, as given 
by I'clas, is identical v.lth that of Greenberg's. He 
splits tt^c orii.'inri LPP in two. T-he first one actually 
optami2.es the cost, wx.ile tt:e second one chooses an 
optimum 0-1 varirble (.a. secuence betv.een h-vo jobs on a 
maebine) out of a set of unassigned binary variables 
(i.e., unresolved sequences). Txiis variable is assigned a 
value 0 or 1 and inserted in tire first part 'which is then 
optimised. Belas hv'=.s simaalified this procedure for solving 
the roalie span problem. He starts 'ivith s. heuristically 
generated schedule. He chooses to reverse the sequence 
between two jobs on a machine which lie adjacent on the 
longest path of the schedule. The resulting sciiedule is 
tested to find vxhetiier we can say that no furtlier decrease 
in the make span time v.ill take place by the reversal :of 
any new sequence in tne schedule. If the tests assure 
it, then a backtracking step is taken and a new branch 
is explored corresponding to reversal of another sequence 
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c longest peth of triat schedule. A backtracking 
tl'.a t brings the scheduler to the root of the generated, 
assures the optinality of the solution from vh.ere the 
nt 1 ; r ek tr Pckijj. r; step vvos originated. This is o, post— 


or,tirait£ tion process. Pll depends upon getting a suitable 
heuristi rally generated initial schedule. 


Pj 1 ihol son (IA) has developed a near optinial 


pi'ocedure to solve the job-shop problem v^ith due dates and 
priorities, his procedure employs a. tliinking which is very 
similar to that of Salas. He also prepares an initial 


schedule 


and then reverses the sequence betv^een the opera 


tiona one 
The rub-0 O' tlH; 
reversals. H 


by one until a lowest value of cost iS' found, 
f'lity is introduced by limiting the sequence — 
e has reported considerable improvement over ta 


■'uristie •‘ik 


thuvib rules. 


Pll the above analytical techniques a-re still 
bagged in the elementary problems. The future of branch 
and bound procedure is very good, because as a.gaanst 
Ba'la 3 ‘siriethocl it gives a procedure to construct tire 
schedule. 
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BIiAl'iCil i'JjD BOOaD METHOD 

The teci-riiauG to solve the integer linear progh** ’ 
ning problerj vos originelly put forvi/ard by Land and 
Roy, , Lomniclii and f e^iv otliers have used it to solve tlic ’ 
shop schoauling with t^jo or three machines. The most 
genersl rornnlc tion of the problem with m jobs and n 
mrchinos hrs been given by Gre-jnberg. This formulation 
is discussed brielly in the first two articles of this 
chapter. The third article proves that making a Gnatt- 
chart is equivalent to solving & linear programming probl- 
(LPP) used by Greenberg to finci the schedule at any stag ; 
of the branch and bound procedure. The fourth article 
enumerctes some of the limitations of Greenberg's method. 
In the fiftn article an alternate procedure is developec, 
on tl^e lines of th^; Greenoerg's formulations. Ttie Irtar 
article gives the bounding proceaure for the uevelopeu, 
metl-oci, 

3 . 1 PROBLM STATMEdT 

Given m jobs and n different kinas of machines, 
each job is independent and. is required to be processed or 
n machines in a fixed sequence, a schedule is to be found 
which is optimal for a given linear objective function. 

The restriction of n processes on each job can be 
dropped, but for simplicity in presentation of the main 



Til ' iine.-^r objective functions can be enumerated 


as follows ; 

(o) Minimising the time i completion of all the 
jobs, roferx’ed to as minimising the make span. 

(b) Minimising the idle time on the machines. 

(c) Given costs of job idle times and machine idle times 
vvhich vary linearly vvith the idle times, minimise the 
total costs of tne schedule. 

’T}.'-e nctitioi.,s used are the s£me as given by Greenberg. 

Let tl'-c' secuence of ms chines for job i be (i=l, 2, . . .m) and 

tl'ic kth el'-menit of this vector be l(i,3i). Let the start 
tiaxes demote the start time of job i on the machine k. 

The a.orfitions denote the processing time of job i on 

hr cl' in e k. 

Then ti'io constraints of sequence of processes on 
any job cm bo \aritten -r.s follows ; 


Sil(i,k)4- (i^k) ^ ^il(i, k 1) 


( 1 ) 


(k = 1.2... M-i) for all i. 

Since only one riiachine of each kino, is available, the 
jobs vvill have to comp-ste for getting cO time slot on 
any machine. 


Let tv-.o jobs i and o competing for a time slot 
on machine k, tlien if job j preccces (denoted by < ) 
job i then, 

Ik " Jit jii 


( 2 ) 



til en 


rr.: if ,1ob i f job 

( 3 ) 

Only one o.f the t.>.o equ-'tions (2) oncl (3) can nold good 
at at tirae for any feasible solution. 


The objective functions, enumera.ted earlier, can 
oe forr-ul "'t'-'cl as following : 

(e) fSiniraising ibc GO’ke span : Add. the following 
constraint 


Sll(i,n)-4-.Di2(l,n) a Sq for ail (4) 

whore Sq is the make span time^ The objective 
function is Min 2 = Sq 

(b) Minimising the idle time on machinefs : Add. the 
constraints of the following type ; 


Sir l-S. 


for all i and k 


where is the tirae of completion of all processes 
on machine k. The objective function is 2^ = min Sj^. 
(c) be can .further extend the objective function by 


a.dding the following constraints : 


"^il (] 


^il(l,n) 


SMjj. -i- for all i and k (7) 

where SJp, SMj^ denote the time of finishing the last 
operation respectively on job i and mo.chine k. If CJj_ 
anil denote cost of per unit idle time respectively 

on job i and machine k, then the objective function can 



bo v;ri tten (io\\n ‘-s lollcws : 


m n m 

L?ir. CcGt = ^ CJi ^ SJ. - b,. F CJi'- 

i=l i=l ^ ^ ^ 


n 

/„ , CM-ib- SMj, -71 >_ CMF D . , 

k=:i ^ k=:l i=l • ^ 


oi , .bin b = .'o Cjj_. SJ. 'i / 
i k 


S^k 


( 8 ) 


Tbo coristrrinta (l) f nd (4) ? or (5), (6) r,nd (7) vvith 

rppropri,f.'te obFsctive function form a linear programming 
p rcbl em . 


24 


1 . 2 unERNBEiKG’S PQRKijL.ATlOI^ Udli^G THE BfiiJMCri MD BOOl^D 
TECbBIfUE 

Tne equ' tions (2) ;.na (3) in 4-1 C're combined 
by Greenberg ac follovvs : 


4,k ^i 3 ' jj 


^ (!-]{?. )M 
-** (3 


where . : 1 if job i 4 job j 

- 0 if job j <, job i 

M = arbitrarily large number. 

If job i < job j then equation ( 9 ; become: 

Q n < 2 ; 

'■'jk 

If job i 3 job j then. 


Sik ^jk 


M 


(9) 


( 10 ) 


Since M is an arbitrary large number, both equations (2) 


and ( 3 ) cm be deduced from equation (9). Addition of 
equation (9) to the LPP mentioned in previous article. 
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mokes it integer LP'f- Let 

31 I 4 - 

-a - x,,)D, 


V ~ ^ ^ Q . . ' Y '''• 

“ '^ik 


ji^ 13c 
.k 


Sub,!<^ct to intep:er restrictions on X.., V cannot have an 

31 ' 

rrhji.rory vslne. In other vi-ords for V to be non-restrictive 

there, shjoiid be no integer restrictions on X^ . A non- 

■^i 

restrictive volne of V implies no restriction is being imp^&r',>d 

on S- anci 0 -, whic-h is equivalent to saying that there are 
'*■ 3£ ■■ 

.more than one machine of kind k. If this be the case the 
constroints (2) and (3) need not be specified. Specifying 


ecu 


ation (2) .'T (3) puts restriction on Sj_}f and Sj]^ and 


thereby pirr.s integer restrictions on X., . Hence integeri- 

31 

nation of X-vari al^le is equivalent to adding to the LPP 
isenti oned in previous article tlie restrictions of type ( 2 ) 
or (3). S-iermfo re Vve need solve a non-integer LPP only. 


There 're (nJ)'^^ sudi X-variables vvhich have to be 

assigned a 0 or 1. The constraints (l) an^L ( 4 ), (5), ( 6 ) 

or (7) being linear and. convex the region '"'f .solutions . c 

is a convex. pbLyhociron in a mui tidimensional space. 

Therefore, tice branch and bound procedure can be applied 

for obtainin,g an optimal assignment of X- variables 

(see the reference 11). This procedure tells that 

the rainimum v;lue of tlie coot function i, obtained by 

I n 

assigning a X-variable a C or 1, is the louer bound for 
£-ny further assignment of yet unassigned X- variables. There- 
fore an assignment tree can be constmeted. based on branching 
from the lower bounds onvvarcis. The nodes in this tree contain 
the value of the cost function 2 . for the pc.rtira assignment 
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.'rici tte br.'’nd''cs ^..-present tiic i-ssignment of an 
X-vrricbj.o or eouivrlently the adaitlon of the constraints 
of typo (2) or (3). Those nodes vvhieh have not been chosen 
for dovinAord. branching are unraarhea. The node with the 
lowest Z a:r!:!onr h-e unmarked node is chosen for downv/ard 
branching, in the rroe. Some branches rc.ay lead to a node 
v.ith no-f ea.sibl e solution. Such nodes ere never considered 
for- c:o'wn'v./rd branching. These nodes indicate that a, 
loopine sequence ha-s been tried upon. 

Proceecing in the above manner will lead to a node 
where al.! the X-varia.blf;S Vvould have been assigned and it 
has the lea.st val.uc of Z in comparison to all other unmarked 
nofi.-s in the assignment tree. This is. tiie optimal solution. 

Selection procedure of node for dovunward branching 
may choose a node higiier up the tree upsetting a. number of 
assignments. Inis is the problem of backtrack and is a 
very serious probleai. This has been considered in article 
3 . 4 . In tiif' next articles ‘‘e will look at the solution of the 
linear programming problem for calculating the value of 
any node in the assignment tree. 

3,3 THE LINEAR PRDGfiijiMliMG PROBLEM 

Greenberg has used simplex metnod to solve the LPP. 
However because of inherent simplicities in the constraints 
the method of obtaining an optimal solution turns out to be 
a simple procedure called the active feasible scheduling 
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procodure. ‘Tlo.is procec^ure nas been used in all the 
heuristic methods of solving the job- shop problem used 
hither to. The schedule tiius obtained is called a Gnstt 
Chart. In ‘this article an equivalence uill be establi- 
shed betv^een the active scheduling procedure and the Dual 
Cl rupl en p VO c ed.u r e . 


fny complete solution to a job- shop problem is 
ch;;- rsctcrissd by a partieulcr assignment of integer 
vr'lued j-varicibleS;, or in other words a particulcir 
sequence of jobs on all machines. Therefore an active 
feasible schedule also has such a particular sequence. 
Subject to such a sequence, an active feasible schedule 
j-h; CO n s t ru c t e d s f o Hows : 

A job is assigned to a machine according to the 
specified, sequences, at the earliest time at Vvhich both 
the machine and the job are free. 

Such B procedure vi/ill generate a schedule in 
which no iricchine is idle for ci, length of time sufficient 
to process completely a scheduls.ble job. 


Let us now examine the LPP 
corresponding to tlie root node of 

HI , n 

Min Z CJ SJ j_ 2. 

i=l k-i 


. Rewriting the LPP 
the assignment tree. 



o r Max Z = -X. CJj_r SfAj. 

i k 
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subject to, S.-, 

1 ) - 


0 


SJ, ^ 


^ ^u(i,n) ^il(i,n) ^ 0 

- Si3^ . ^ ^ 


( 11 ) 


0 for all i ^ 

The coefficients rr+T.- ^ 

- -iix of the LPP coht-ir. 1 
sither 1, 0 or -i rr, . ^"-ins elements, 

‘*'^ere are a large m . 

tt-is ■r.atriy. of i-eros in 

,, ^"'^'^ipulrtion of ev- ■ 

the esGenti-l esrf . f . . ^ this matrix is 

‘-■’1 simplex or Duel .-• 

This is enorTnu'-' -r. ^*^plex. procedure. 

" ■■ort oven for a sm-h-. 

(11) tile si^.e of tr^- problem. In 

coeliicients matrix ; 

X i'-pq vull be like. 


P (h-1 ) ^ 

^ - X mn 


■'Tid 


0 = in r m . . 


T*:.. r.u..pl,3 .na Integen-f 

Of tyne (2) or C'-o , ^r^o.t,ion eouations 

l-') v-m iperease this si.„ 

^^4,6 lurtlier. 

P'urther ev 

-hetion of the cons + >,, • 

rs given in (11) , temts of the LPP 

’ revecls that starting ti 

an t Job is aepeMpn* ., operrtio 

^tnr onpy upon the r-o - 

epor.' tion ftliich clirrr.+-i * r>ish time of the 

- oirectly precedes it. t i. 
constraints of tyoe i'ss . . ^ reduction of 

■ U) or (3) introduces tb„ s . 

precedence relation c. inter job 

on B machine. Hn' 

is sn obvious flovk lit ^ s,.ows that tliere 

i-iKe structure of tf 

sciieduling procedure oxp^oit^ thi Active 

xp^oito this property. 
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In the Du Pi simplex procedure for a gener&l form 

of LLP^ 

Mcix 1 r c X 

subiect to A X >■ b 


the necesssry end sufficient conditions for obtaining an 
optimcl solution are as follows : 

X for basic variables (12) 


t-, •it’y 


- ^-fii variables 


whore C_ 3 ^ : cost of ith basic variables 


( 13 ) 


C-r 


cost of j th non-basic or basic variable 


and. (Y^p : the coeflicient matrix in the tables at any 
i teration. 

In the LPP for tlie iob-shop problem all these costs are 
negative. 


To esttblish tne equivalence in the Dual simplex 
and the,' active scheduling procedure, following theorem can 
be stated. 

3liehr^_l: F'or the given set of constraints the active 
►-> cn ecLui. mg is cn optiiui'.'l schedule, i.e. it satisfies 
equations (12) and (13). 

Proof 

In rn active schedule a job is assigned to a machine 
at the earliest time at which both machine. cind.'.j-ob- are, free. 
Therefore all the variables in the constraints are greater 
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or equol tc zero. So the active schedule is a feasible 
schedule. This satisfies the criterion (12). 

The price varifbles SJ^ (i = l,2...rQ) and SMj^ (l!:=:l,2. .n) 
and all the start- time variables except a few have 

positive nonzero values so they \fkill be in the set of basic 
variables. The nonzero time gaps between tjiy two operations on 
a job or a machine will also be included in the basic set. 

These are tiie surplus variables in the Dual simplex fomulaticn 

fny row in the optimal tableu can be written as follows 

Si4- Z b = 

all j 


vvhere 

Sj_ ; basic variables 

: value of the basic variables 


(Xp^^) _ ; those for any i or surplus 

variables which have zero values, i.e. , 
the nonbasic variables. 

For the active schedules we can also v/rite similar equations 
for the nonzero variables in terms of zero valued 

I 

variables (Xp^) ^ . 


^i 


y Y. (x ) = b^ 

^ pq j ^ 

all j 


(14) 


L'l' is will be a tableu similar to the form obtained in th( 




Increesing (^pq) ^ above i.ero level is equivalent to 
pu siiix'ig a. job to the right of ti:e mlninam time when both 
tx- e mach-ne and job arc froe^ Such an assignment can push 
some or ell the basic variables to have more positive 
vaiut. th.rii «vhat they uhll have if all the assignments are 


done at the earlipst possible times. This 

a situation where in equation (14) the Y. . 

^ j 

0 or -1. 


corresponds to 
are either 


Then tile optimality criterion (13) ? 

^■'j ~ ^'ij ~ nonbasic j 


becciusG all the costs aire negative and Ypj are either 
0 or -1. vSihee both the criteria (12) and ( 13 ) are 
satisfied by tiie active schedule;, it is an optimal 
schedule, ilite teat it is optimal with respect to the 


corstr-'ints given in the corresponding LPP. 


QED. 

Addition of a new integer! zing equation of the 
type (2) or (3) to a LPP may make the present solution 
to be non feasible. The additional equation can b-e 
auc'iod to tine tableP of the dual simplex procedure after 
exprc'sslng rll the basic variables in it by the nonbasic 
variables of the original LPP and including the surplus 
variable into tine basic set. 

If the additional equation is 

^ik 


(14) 



introducing the surplus vni^lable x, it becomes 


lI 




X = D 


3k 


( 14 ) 


If S-;v f'bd S^v are the hasic variables in the 

-Lit J 

original LPP^ then using the equations similar to (13), 
they can be expressed in terms of the nonbasic va.riebles. 
Then surplus variable x bill be included in the basic set. 
Since nonbasic variables do not have any price the 2^-Cj 0 
r cm. a ins intact. 


If job i and job j are not in clash over the 
laachine k as doterrained by the solution to the original 
LPP, tiien x bill have positive value otherwise negative. 
Hence the solution to the original LPP togetlier with the 
neb ba.sic variable x is a non-feasible vuith all ^ 0. 

Therefore it can be used as initial solution for the 
iterations on the nev».LPP. The post optimisation over tlii 
initial solution v\(ill require & few iterations. 

The active schedule for the new LPP will be 
constructed as follows : 

(l) Add the adaitional constraint on sequence to. the 
assignment tree. Mark all the operations as 
'’unassigried'h 

(.2) Search the tree to find a leftmost operation on a 
job which proceeds all otlier remaining unassigned 
operations on all the jobs. Assign it to the 
appropriate machine a^t the earliest time a^-t which 
both the machine and the job are free. 



(3) Ikev'k tills operation as ’’assigned”. If all operations 
are assigned the schedule is complete^ othervvise go to 
'tlee step '(2). 

This procedure vvill cost more computational time 
than tile post optimization in the dual simplex procedure. 
liov.'Gver in the brainch and bound procedure any backtracldng 
Vvill require construction of the nevi.' tableu by the search 
of the appropriate portion of the assignment tree because 
the storing of tableu corresponding to each' node is a huge 
effort. It Vvill require considercible manipulation tiian 
fo rasing an active schedule. Mother point in favour of the 
?ctiVf.‘ schedule is tiic.t it need store only the assignment 
tree vliile the -dual Siraplex procedure vvill have to store 
tho coefficient matrix or the Tabloj in addition to the 
a s fj i gnm en t tree. 

Storing a Icirge matrix is a difficult problem on 
the presently available computers. Osing backup stores 
like disks or tipes Vvill involve mechanical iroveraents for 
any manipulation in the matrix vvhich vvill be very costly. 
Therefore the choice rests v.itii the active scheduling 
procedu re. 

A further reduction in the calculations of active 
scheduling, i.e. preparing Gnatt chart, Is possible by 
modifying the brancli and bound procedure used by Greenberg. 
It V'.ill be dealt in a later article. o 
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3 . 4 r.T^^l TATIOrjS OF TriE GKEEI\lijEfiG ’ S METHO D 

i-vt fiiy stigfc: ol scheauling, Greenberg fixes en 
integer vcrirbl« to be eero or one, vvhich is equivalent 
to fixing the rel r.tive sequence between tvvo jobs on a 
nr chine. If there are m iobs, then each job requires 


m r c' 

hine. I 

f tl' 

■.lerc c. 

fiz 

‘-'tion of 

Co- 

-1) se 

S Cl- 

euule is 

Til 

K.ed - 

ope 

ration c 

an 1 

DC 

the 

1 operr-ti 

ors 

preca 

\i,e 

have.' see 

:n i: 

n tl.'.e 

thr. 

' nedive 

S c'i"-. 

edu 1 e 

ti.'.c 

;> tree ed 

'Ch 

time ■ 

So 

the constru 

ction 


lest article in the construction of 


time 


con SUIT' in g. 

In the branch and bound procedure one has to bacK 
tracK to all those 'unmarted' (or unbrrnched) nodes *1* 
have a cost lower than dost of the optimal schedule- 

li,e hach-tracMng ear. be minimised by fixing those seque^ees 
earlier which give large difference between the lower m 
the upper bounds. Greenberg fixes the sequences in a 
completely arbitrary menner, end also tnose sequences ^ 
fixed in which the jobs concerned have no clash. BesW®" t 

,, c.-.-n hp pxDected to cause only 
fixing of one sequence cen be cxpecee 

■ cost or the objective function, 

small increase in tne cost qi. ta 

bP-hvwP(^n tbe lovter and tiie upper 
Also the difference beti^vecn xx*c j-u 
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boimes v,ill xiot be sul fici ^ ^ 

to stop the backtracking. 

“■lor-eover 'we cio not expect x, ^ 

^--t a 30b Vvill be sequenced 

earlier to a job which isri-^o.. x, n- xx 

^■^'■cea much earlier than it 

in the schedule. Tnereforp fs ■ ■, 

iixing such a sequences is 

just a futile effort. So tbp r^ x.-, n x i • 

problem of be.cktracking 

is really big. 

-Ano trier limitatiori nf i ., ■, x- ^ t • 

- ui breenberg’s metiiod lies 

in the calcul^tion of the pbo-i- ^ xi x- t x i t 

^ costs of the partial schedule, 

or the bounds. i'ccorciinp tn u-,- x t 

s 10 his metliod the bounds are 


calculated from tl'.e follovi/ 


equation 


m 


iml 


n 


CJ^ Sojh- Q 

k=l ^ 


(15) 


vl.ore SJ'i = time of finishing the last operation on 

job p. and oivijj. - time of finishing the last operaion on 

machine k. These times are based on the schedule with 

partial fixation of sequences, if ^11 the jobs on a 

machine arc in clf.sh witl- pp^v- xx xi xx n-n u-iTt 

fcach other than the Sh]^ vail 

be oven less tiian the sum or xu • „ -ti 

^ -►i fPe processing times of c.hi 

jobs, on tiia.t machine. 


P.11 these Iteitctions tell us that several 
Tsodificatlons ire needed Greenberg’s branch and 

bound procedure. The todiacstlons are suggested in 
the following articles serially 
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3-5 I'-IXI^IG OF :.:iORE OKE SEgOEl^CE AT A TIME 

It is necessary to fix more ttian one sequence at 
8 ti3Tie tile lo-v^er bound colculated should be nearer to the 
probable optimal cost, to reduce the backtraclring. Tvvo 
methods can be suggested cs folloiivs : 

1) lUx the sequence irl" en unassigned job vdth 
respect to :all the otlier unassigned jobs r^r-uiring 
processing on tlie same machine. If tiiere are r such jobs 
then r nodes must be celculated, each node specifying the 
job tiiat proceeds all others. Since one of these 'v.ill 
definitely be assigned the first place among these jobs 
in the prob?.ble optimal schedule, therefore, it viill give 
us tl.' o 1 o i(v e r b o u n d , 

It is equivalent to assigning a value 1 to r integer 
variables as against to only one in Greenberg's method. It 
vtoulc! have taken calculation of 2r nodes plus the necessary 
backtracking to fix c'n operc-tion like the above. The price 
is paid in terms of tne fact that if backtracking becomes 
necessary a greater amount of calculation is needed. In 
greenberg’s proceaure backtracking to an earlier stage 
changes only a few sequence, while in tlriis method all the 
r sequence, while in this method all the r sequences are 
upset. This price is not much because the difference 
between tlie lower bound and the upper bounds in this method 
is larger than those in Greenberg's method. 
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Ii-ti iH'-jj r V i s ol t^ii s niGthod. ov©r "th© Gr©snb©rg’s 
procedure v,ill be further eleborc.ted after presenting the 
s e ccmd no ci i f i c a t io n . 

f.) Consicer the- irQcedi&tely next unassigned operations 
on all j')os. C^'lculate t;..© nodes correspondings to fixing 
each such operation er'rlier to any other opera..tion requiring 
srv-e machine. If there are m such jobs then m nodes vvill 
be calculatoc. One of these operations vvill definitely be 
placed f.t its decicra position, therefore it will give tlie 
luWFi- I'-jund of tl.'o scnedule. 

If i'h tic-n on a machine k is chosen to be 

assigneci prior t .■ other(r-l) unassigned operctions on the_ 
pane maiohine, then the above process is equivalent to 
fixing r integer values to one. Hence as regards its 
comparison vdth Greenberg ‘s method the arguments given in 
tiae first meti-od. also hold here. 

i'S we can see th.ere is an automatic selection of 
next cl ashes to he resolvee in the above metiiod. In tlie 
first method there is a problem of finding a suitable 
machino on Vvhich. the clashes must be resolved. This con 
be solved by choosing a raachine whicl'’. requires processing 
by atleast one of ttae iKimedlatoly next unassigned operations 
on all thejobs. Such a set of operations can be called a 
MIX. 
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vvi liii iTai-i-s V. (.■) will co!Tipc''.rG ths t’l/io TistHod.s by th.© 
help Df en oxrmplG. Consider r following problem 


Job 

1 

1(4) 

2(3) 

The costs of job lateness 

Job 

2 

2(3) 

1(2) 

and machine idle time are 

Job 


1(1) 

-f2) 

unity. 


The numbers inside brsckets are the processing 


timos .'•nci outside the brockets are the machine number. 
ii--a i.natt on art in ti'-e first procedure is os iollows : 


jj 

--jr. 

X. li., ; ~ j 

' 1, , * 




1 ; h . . 




Oj 


f 4 -) 


/Ju) ' 



^ -7 1 

._.xi 

3 1 

I" 

V. i 



, 

1 1 ■ *! * 

r n /. 



■i ‘ 7i 


hh) 


1 ^ 

iX ' -rp ai 


1 " 3 i 

. ): vl 


:'2-' (' ^ ;t! 




'! ) 


'€^'' 777 ^ 7 :. 

I •> ; 


1 


TiiO superscript is tlie machine number. 

The subscript is the actual time of the 
completing tiie operation. 


Fig.l 
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.jk' : Operi'tior! ik directly preceeds jk. 

For tl.'c C' leu lot, ion of the root node;, see the Gnatt 
chart No. (1) in the Fig.l. 


SM^ =7 SJi = 7 

SM 2 _ = £ SJ 2 r 5 

SJ3 = 3 

Liachine 1 : It is obvious that tiie operation 31 should 
not be placed as the first job on machine 1 
sequences : 11<”< 31 21 ; ^ - 6 

.11 /C 21 ./< 31 • A 1 - 7 

v- t..^2 ~ 

31 </ 11 3' 21 , - 3 

31 C 21 /- 11 ; /'i _ 5 

4 

= min (Ai 3 ^ -'4) = 3 

Machine 2 : It is obvious that operation 1.2 should not 
be placed as the first job on the machine 2 


sequences : 22 32 12 

II 

32 A; 22 3' 12 

; . 5 

22 <^<12 <^, 32 

; = 6 

32 <(<12 <A 22 

; O4 = ^ 

= 3 


Since v/e have not considered how does a sequence on one 


machine affects the sequence on the other machine, we 
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ccn'^iot PC... ti:c Pnu however, tiie largest of 

cTiCi /\^' nill be the r.inimijra increase in the cost function Z. 
Hence ^ - r!ax( _^^ , a-^) - 3 

4 = S% -- SJ 3 _ 4-- SJ 2 4- SJ 3 


y\ 


- '/ -£ 


g -i- 7 -+- 5 +-3 +3 33 . 


Tl:-o fix in Gnatt chart (1) has tire operations 11, 22 and 31 
.60 tl., nPciPines 1 and 2 rre involved in this mix. Ihe 
.jt';iOlioci 3I p-i-Ves a pusi'i oi 1 unit to other operations on 
t'.c no: .c nieu .ai.so tJce least push given by any 

oi'ior operatio.rj on vaci'-ine 1, in the ruix. The minimum push 
or: :r-,-c,hinc 2, -as given by operation 22 is 2 units. Let us 
then choose the machine 2 for resolving tiie clash for the 
first plrcin.r: on it, becouse by a larger push it is more 
1 lively to ease clasi-es among other operations. 

On the machine 2 only jobs 2 and 3 are thie 
contenders I'or the first place. 


tion 

22 in 

firs 

t place 


See 

t.£ e ( 

hiatt 

chart No. 

(2). 


mi 

= 7 

SJ^ 

.. 7 


foidl 2 

r 8 

SJj 

= 5 




SJ, 

= 5 


Machine 1 : 11 yh 31 <(4 21 ; — 1 = 

11 « 21 4 31 I l\^ = 

31 « 11 21 i ^ 


31 < < 21 24 11 j A ^ = 


6 

7 

3 

5 
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n pPTv^t-.inn ^<2 in first vlc-ce 

See tile Gjirtt chert ao.(3). 

r 8 8J3_ = V SJ^ = 8 

Sih; =8 SJo = 3 

l'/> eclaric 1 : ll4<^31C-4 21 } A 2 = 4 



31.^114421 ; 

1 

= 1 



ii 

Machine 2 : 

32 // 22 // 12 ; 

\> 

It 


32 « 12 </ 22 ; 

n" 3 4 


/\^ = .2 
A = 2 


4 r 36 

The operation 22 has tl'ie lower hound so it will be placed 
as the first joh on machine 2. 

The new MIX : 11, 21 and 31* 

Machine Chosen : Machine 1 (only one machine in the mix) 
The operation contending fo^^ first place on machine i 
are 11 and 31. 
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0; ..'P.' i .101 


4 4 .-iv. T ‘i 


' O}’- 


411ilLJoib_on r.^coiiiio 1 


.?5 ;■' chine 


js. 1 X . 


Si the Gi?.- tt ci'-crt Ko. (4) 
hhi, = 7 SJ-j_ = 7 

■^'4 = ^ 

1 1 1.1 .V. II 31 

11 "h 31. ■• 21 


I'p- '7- ''<32 

i-.' </ :<„(/ !•' 


J) - 




SJ2 




= 

-^2 

= 


1 

1 

2 

3 


A= 2 


Z = 38. 


A 2 
/ \ 


6 


‘iiL. .‘ILl'i 2.!l...._3l o.'’ ti- !.:i l'i r:.h.. .lob on me. chine 1 

SU e tJ-e Gns.t!'. cx.:Ci-t ]:^o.(5). 

®:3_ = 7 SJ^ =8 SJ2 = 5 

Sj72 ^ S SJ 3 = 5 

.?Af chine 1 : 31 « 11 21 , = 3 

31 << 21 1.L ; ,..4\ = 4 

a1 

A = 2 

chine 2 : 22 A.32 a, 12 i' /\^ = 0 

4^1 ~ Z , Z « 3 5 • 

The o’perptlon 31 has tee lovser bound so it v-vill be placed 


£s tile .first job on machine 1 



In tie Gnat^ cnert corresponcing to operation 31 
v.e observe that 32 is pu tomatically scheaule as it is 
the ney. OKer. tion alter 31, it has no clash v.itii the 
rs’iifining operation, on the second machine, i.e. , the 
or or.: tion 12. 

The noH nix ; 11 anh 21 

Th e TE a ch- i n o. cl lo s on : 1,^ a ch i n e 1 

The contenccrs for 2ncj, place on Machine 1 are 11 and 21. 


Operation 11 in the s econa piec e 
See tile unatt cl'X'rt i'to.(6). 

^ = 7 bJj^ - 8 SJ2 = 7 

GL^., - c 13=5 

a r cl : i ( i 1 : Id,, c 1 '• s I‘i fi s , - 0 

Machine 2 ; ho clashes, ilf' = 0 

1 = 35 

Operation 21 in the second place 

See the Gnatt chart ho. (7). 

= 9 SJ^ - 12 SJ2 = 5 
® f>M2 SJ3 = 5 

Maciiine 1 and 2 : I'Jo clashes 

h 3 43 

The oporotion 11 is schecoleu in the second piece on the 
me chine 1. i/uitii tliis all the clashes are resolved so the 
scheduling stops, hote thr t 4 is not tl'ie nctuel cost. 

The cost is calculated as the sura of the job-idle times 
and machine iale- times multiplied by their cost per unit 
time. The Optimal cost is 5- The assignment tree is 


given in Fig. 2, 
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Tbis I'Tobles: r.es been solved by H. Greenberg (10) to 
miniraise only the idle- tines on the machines in 19 nodes 
Tif ter solving 19 linear prograinining problems. 

Ti.e cclculotion of bounds for c general mxn shop 
problem '-re discussed in greater detail in article 6 of 
this chapter. 

Fortunatelyj, no baclitraclEing is required in this 
problem, if backtraclcing becomes necessary then a search 
of the tree becomes necessary to find the node from where 
the further branching vdll be resumed, ihe Gnatt chart 
will be constructed according to the sequences fixed Vvhile 
traversing the path from tne root node to the chosen node. 


The assignment tree for the same problem using the 
second method is given in Fig. 3. Tne bounds are calculated 
in the same way as ticat in the first method. 




4-6 


The dotted po/tli shovvs tlie branch and bound path 
adopted by ti:e scheduler. In this method we calculate 
three nodes at each stage of the tree while, in the first 
method this is not & fixed number; it can be at most 3 
equal to the number of jobs. Note that in the above tree 
nodes A cUcg B are identical yet they require to be dupli- " 
cated. This situation arises because at the first stage 
of tree we have compared nodes marked 22 end 31, i.e., 
those operations which require operation on aifferent 
irischines. There is a. aefinite possibility of both being 
s eh ecu led to goth. or. 

This prompts us to think of scheauling more than 
one operation out of the mix of the imroediately next un- 
asslgned operations on all jobs, the one we talked about 
earlier. This process will be equivalent to simultane- 
ously attending the queues on different machines Vvhich 
occur in the mix. It can be easily argued that we can 
be completely sure that fixing one opera.tion in the 

mix will give e lower bound but not so about fixing more 
than one operaition in the mix even though they require 
processing on different machines. For that then we hove 
to consider the combinations a-mongst all those operations, 
not necessarily in the mix, on various machines thc».t are 
in the mix, ond are in clash, with each other. This 
may swell to a large number. The number of nodes 



corresponding to c.ll the co.-sbinations of tnem vdll be 
even larger. So it is better to leave it. 

Let us come bock to the comparison between the 
first and the- second method. It is quite possible that 
the operation corresponding to the lowest bound in the 
first method may turn out to be one whose predecessor 
operation is still in clash with other operations, that 
is, it is not in tiie raix. If the predecessor is shifted 
to the right as a result of some, clash resolution then 
the earlier calculation may be upset, i.e., the back- 
tracking may become necessary. Besiues this, the Gnatt 
chert construction becomes difficult. In the second 
method each node can be labelled with the starting time 
of the operation concerned. In the first method we Vvill 
have to search the tree to find out how many operations 
to tile right of tile recently sdieduled operations get 
pushed to the right. However, there is a way out if we 
utilise the idoo o'f a Mix, If any operation corresponding 
to the lowest bound is not •' member of the Mix, then 
keep this modiine aside. Choose another machine from the 
laix and again find the lovivest node and repeat the above 
test. We will definitely find one operation, as it is 
assured in the second method. Ihis will save us for 
complex Gnatt cJiart construction for each node. 1316 price 
is paiid in terii.s of the calculation of a larger number of 
nodes. However, it maiy happen that we do not have to 



consider eny ne'Cs node. If ell the operations in a mix 
require operstion on a single machine then no extra 
calculations are needed. If there are more than one 
machine in tiie mix, then there is a possibility of 
simultaneous scheduling of one operation on each of 
these machines. Therefore, it is a safe r bet against 
large computer time in constructing Gnatt charts for all 
the nodes in the tree. Moreover, vvith tile new arrangement 
in tile first method, v.e 'will be able to select 
a clash which gives large shifting of the scheaule to the 
rignt. This is likely to ease clashes among the remaining 
oper' tions. This is not possible in the second method, 
because the choice of the lowest bound may fall upon an 
operation which gives least increase in the cost. The 
least increase in the cost will usually produce small 
right shifts in the schedule. . 

From tile point of vievc of the number of nodes to 
be calculated at any stage in the tree and the duplication 
of nodes, the first metnod weighs higher tiian the second 
one. So from now on we. will restrict our attention to the 
first metiiod- For ease of reference we can refer it ais 
the alternate procedure in furtiier discussion. 

In the Gnatt chart no. (l) of Fig.l we have seen 
tiiat for deciding the first place on machine 2 we 
cor-usidered only jobs 2 and 3. Out of all the three jobs. 
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tne earliest time of finishing is for the job 3. Both 
jobs 2 and 3 have' the start timesless thou this time. The 
operation 12, on job 1 has a starting time of k, vh.ich is 
larger than the earliest finishing time. So we Isnow i^at 
as far as .the first placing on machine 2 is- concerned, 
the job 1 is not a contender. This observation can b4 
generalised as follows : 

for deciding sny placing on a machine, only those 
unassigned operations aiust be considered Vkhich have their 
stert*^ time less than the time of finishing an unassigned 
operation which has the least time of finishing in tlie 
present 'Gnatt chart. 

This Viill reduce tii'- number of node calculations 
and the si2,e of the assignment tree. Therefore, the number 
of the bccl<’tra.c>:ing sttvps can cdso be expected to decrease. 

'vve are now in s position to put dov.-n the alternate 
procedure in tloe form of an algorithm as follows : 

(1) Form the mix of ''leftmost unassigned opera? tions on 
all jobs. If any operation in the mix is not in 
clash, then declare it 'assigned ' and replace its ■ 
position in the mix by its successor operation. 

(2) •Choose a machine out of this mix, operations on vvhlch 
produce the largest of the minimum right shifts in 
otlfier operations in the Gnatt chart. 

(3}- Calculate the bounds corresponding to those operations 
on this machine v^hich are in ' clash' ^ as discussed 
earlier. If the lowest of these bound is corresponding 
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to pji operation Viihlch is not In the then choose 

another nachlne according to step ( 2 ) and calculate the 
neVi' nodes. Repeat this process till a lowest bound is 
found corresponding to the operation in the mix. 

(1) Search the lo Vi,' er bound among tne nevdy calculated nodes 
and all the unbrajiched nodes in the tree. Further 
branching is done from this node onwards. Update the tree. 

A backtracking step at this stage may lead to a node 
the rssociatea operation of which is ahead of its corresponding 
mix. In such a case go to step (3) where the new clashes are 
resolved on a mcchine that requires processing by the operation 
which is preceding tx^e currently assigned operation and is In 
the current mix. In subsequent branchings Gnatt-charts are 
constructed such that no seqiences are violvated tlaat have already 
been fixed in the current branch. This procedure is followed 
only till there remains no assigned operation vliich is ahead of 
mix. 

If the above change is not required go to step ( 1 ). Ihe 
procedure is stopped 'when all the cle^shes are resolved. 

Rext we shall compare the alternate method and 
the modified Greenberg’s procedure. The idee, of mix can 
also be introduced in the Greenberg’s method. The 
modified procedure is as follows : 

(1) From a mix of tiiose leftmost ‘ unassigned operations 
on each 30 b. The unassigned operations are 
those Vi/hose sequences v^ith all otiier operations 
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on the S'CEe ma chine £.,re not resolved .completely. 

If sucii an unas signed, operation does not have any 1 
clashes then it is removed from theimix^and its 
successor is placed in its position. 

(2) Choose 8 mechine in tJr:e raix in the same manner as 
In the alternate procedure. On this machine choose 
tv^n operations vdiich produces the largest minimum 
shifts rela.tive to each other. Calculate the nodes 
corresponding to tifno seauences for these operations 
and calcul^'te the bounds. The sequence with the lower 
bound is selectea. 

(3) Dpdate the mix and the tree. If all clashes in the 
scnedule are not resolved go to step (2), ottlerv-iflse 
stop. 

Let us then solve tine 3x2 problem, solved earlier 
by the alternate procedure. 

The root node is same as that calculated in the 
alternate procedure. 

The macij-ine chosen from the mix is likewise 2. 

Tine operations in clash over this machine are 
22 anc. 32. 

Following ere tl-ec'al-culations of ttiose nodes 
which were not calculated in the alternate procedure. 

Rest of the calculations are left to avoid reijetition. 
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^equ ence 22 32 : See Gnatt chi-'rt Ko. (2) in Fig.(l). 

2 =35 

Sequence 32 <^ 22 : See Gnatt chart ho. (3) in Fig. (1). 

2 r. 36 

The node marlsed 22<^ 32 is selected. 

Kew mix : 11, 21, 31 
Machine Chosen : Machine 1 
Operations Qiosen : ,11 and 31 

Sequence 1 1 <^3 1: See Gnatt chert ho. (L) in Fig. 4 


SMi = 7 SJj^ 

= 

7 SJ 2 - 5 

SM^ = 9 SJ 3 


7 

Machine 1 : 11 ^31 <<.21 

> 

= 2 

11 <^21 << 31 

» 

y 

4.1 

21 <.11 « 31 


^^^3 = 10 



= 1 

Maci-iine 2 : 22 <^<32 < 4 12 

» 

y 

Af = 3 

22 <<12 <<32 

m 

y 

A|= 2 


^ r 37 

Sequence 31 11 : Sec Gnatt cn: rt Ko. (4) in Tig.l. 4= 35. 

The node marked 31 <Cll is selected, 
in tile Gnatt chart corresponding to this node, th 
operations 31, 32 anu 22 have no cla^shes so they are 
declared ''scheduled” as far as this branch of the tree 
is concerned* 

The new mix : 11 and 21. 



Secnence llv 21 * Onc t-t chcvrt Xn Fig.!. 

S = 35. 

Sequonce 21<^-I1 : See Gnatt ciiart Ko. (7) in Fig.l. 



Fig. 4 
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O-’T liV, ,,. 


. V -2! <4 
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Fig. 5 

in considering the sequences on various machines 
for the calculation of we have to he very cautious not 
to violate tiie already resolved sequences on the tree. 

For this we must keep track of the resolved sequences on 
each operation indiviaually. In the alternate procedure 
we need track only the scheduled*' operations which is 
relatively a smaller job. 
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Ihis taetfetod S2.-&0 gives the solutioH iyi sev^n 
It SO happens ‘becc.use the ^obs fit in very viell Vi/ith each 
other in this problem. 

In constructing the Gnatt chart in the modified 
Greenberg’s procedure, one has to be cautious that no 
previo\.isly fixed sequences are violated. Therefore, every 
chart must be constructed from the beginning maintaining 
the sequences already fixed. In the alternate procedure 
we need only reform the previously obtained chart by 
pushing all those unassigned operations which are in clash 
with ti'io selected operation. Only when a backtracfefcng 
step is required the Gnatt chart is built from the 
beginning. Hence the construction of chart in the 
alternate procedure is very easy. 

Let us compare for a general case the number of 
nodes which must be calculated to solve a problem. Suppose 
on certain machine i-n tine mix there are r jobs in clash 
with each other. Note we shall consider only those jobs 
to be in clash vvhich are in clash with an operation willi 
the earliest completition time. The alternate method will 
require, 

r i-(r-l) 4 (r~2) 4- . . . ■4-3-V-2 - ^ - 1 

number Q.f nodes to deterralne the complete sequence of 
r jobs. 



In tlie modified. Greenberg procedure, Tor- the 


iob v.i:lch vd.ll be pieced first, sequences with r-1 ^obs 
must be rer--olvv:d. Out of the remaining r-1 jobs, the 
next ]'}.iecing is decider, in e t cost r-2 sequences. Hovvever 
If the resolving of earlier r-1 sequences remove clash 
among some of these r-1 operations then the sequences to 
be resolved next will be less than (r-2). Ihe maximum 
number of sequences to be resolved is as follows , 

(r-1) (r-2) 4'(r-3) -V-. . . -v*3-t2-tl - 

The corresponding number of nodes is r( r-1). 

The minimum number of nodes can be easily argued and it 
is 2(r-l). The number of actual nodes calculated depends 
on the combinations of ttve processing times, the costs of 
job-lateness and the selection of the operations and the 
mechine for resolving ti've clashes. 

A table for different values of r is given 
for comparison. 



r=2 

II 

r=4 

r*:5 

r=:6 

r(r-^l) 

-1 2 

5 

9 

M 

20 

rXr-l) 

2 

6 

12 

20 

30 

2 (r-1) 

2 

4 

6 

6 

10 


The alternate method nas a number which is the exact 
ari thraatlcal mean of tiie maximum and. the miEimum numb.er 



of in th-o moflifiod Greenberg *s procedure. For 

s’lterne'te procedure tiiis riuniber is independent of tiie 
coiTibinatorial structure of the problem. 

Ihese numbers are little deceptive because they 
strongly depend upon the baclttraclsing. The backtracMng 
in the alternate method vvill upset more sequences then 
in tile modified Greenberg’s procedure. However, because 
of fixing (in the alternate metliod) the sequences amongst 
more than one jobs all together increases the discriminatory 
povver of the bounds and so the number of backtracMng steps 
in this metiiod can be expected to be smaller than in its 
counterpart. 

Ti':.e above comparison establishes that the alternate 
method is relatively less aepndent on tile combinatorial 
nature of ,ti?e problem and the price paid for it is also 
only moderate. For the problems like dob- shop we cannot 
define a closed form formula to prove the superiority of 
one procedure over another, 'we can compare then only in 
terms of the time taken for solving and the amount of 
dependence of tile procedure onto the structure of tiie 
problem data. 

3.6 THE iliTERilATE METHOD AKD THE CRITICAL PATH METHOD 

The scheduling of m jobs on n machines can be 
considered as a project. The critical path method (CPM) 
can be used to construct tiie tree of seouences for this 
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pi>a^Get* This tree is constructed in much the same 
raanner as the alternate method does. The node represents 
the fixation of an operation on some ;3ob. The CPM enumerates 
ell the branches leading from the start to the end. The path 
that gives the least cost is the optimal solution. It is 
essentially a complete enumeration procedure. The alternate 
method discards at every stage of the tree, those nodes which 
have bounds (costs) higher tlaan tne least bound among all the 
unbranched nodes in ttie tree. Thus the siae of the tree 
in the alternate method will atmost be equal to that of the 
CPM but on an average it Vvill be much smaller than the 
later. Increase in the number of jobs increases the size of 
ttie tree in CPM at a much faster rate than in the alternate 
method. Hence tiie •.Iternete method vvill be much faster' 
than the CPM. 

3.P THE B002-SrS 

The bounds are the evaluations of the cost function 
Z resulting from the partial assignment of the sequences. 

The node wi til the least bound in the assignment tree is 
chosen for further branching dowmvard. 

. BrortTi, LomnicM, J. J3.D. Gupte (15) and others while 
worMng on simpler job-shop problans have reported that 
the amount of backtracking depends upon how good the lovver 
bounds are. However, tlie larger the difference between the 
lower and the upper bounds, the lesser is tlrie proba-bility 
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of reversing the choice. Iherefore not only the la^er 
bounds, but also the difference between the bounds is an 
iisportant criteria to reduce the backtracking. 

As mentioned in article 4 of this chapter, ■»€ must 
evaluate the tiie time of completion of the last operation 

on the ktii machine more cautiously. Brown, Lomnicki and 
J.S.D. Gupta have used for the time of finishing the 
last operation on kth machine among the so far scheduled 
operations plus the sum of the processing times of other 
remaining operations on this dadoine. Ihey have not consi- 
dered tiie clustering of these ^obs. For example consider an 
incomplete schedule, in v^ihich the start times and durations of 
unassigned operations on the kth machine are as follows : 

JOB 1 JOB 2 JOB 3 JOB 4 
Starting time 0 .1 8 ‘ 10 

Duration 4 1 3 5 

Brown’s bound will give, 

SMjj. z 4 1 3 5 = 13. 

If we draw a Gnatt chart and add to the above sum the un- 
avoidable slack times, the SMjj comes out to be 16. Ihis 
difference will be oohanced in a problem with large number 
of jobs and non-identical orders of processing of different 
jobs on the various machines. 

Since it is cumbersome to represent the calculation 
of SMjs; by a mathematical expression, it will be presented 
by flowchart as follows : 
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All the lobs rofarped to in the following cfearts 
are the unossigned or pertielly assigned jobs, tha.t is, 
those :obs dbose sequences on a machine are not completely 
determined. 

(l) Calculation of SMjj : 




t 
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A faster chart is as follows : ‘ 




In this flowchart we need not prepare the list according 
to increasing order of start times of the operations ik's. In 
preparing the list Vve will have to search as many minimum start 
times as the number of unassigned operations. In the above 
chart if all trie operations on k are in clash with some one 
or other, tiien we have to find the minimum only' once. 
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Next we will estimate the effect of S.1]^- on the 
schedule. If no oper: tion on machine i: has a finish-time 




•hVi .-r-in 


ell trie unassigned operations must be 


considered, .for possible shifting to the Sl]^ or more 
accurately efch of the unassigned operations must be tried 
as a last ;iob on the kth machine. The resulting minimum 
increase in the cost, can be addea to the cost function Z. 
This is tile ipproaou J.E.D. Gupta and others have adopted 
in their methods to solve a less general scheduling problem. 


Trying ec'ch operation as the last operation on 
mrcioine h af '‘ects trie ^ob completion time of only one oob. 

If there is clustering of unaissigned operations in the 
lel'u end of tlie sciiedulej then Vv'e can get a higher increase 
in the minimum cost by trying each operation to be scheduled 
earlier to others. This v»ill shift more than one job 
towards right and thus affecting schedule in a wider sense. 
Let the minimum increase in cost beA^. 

If an unessigned operation at the left end of tiie 
schedule does not clash witn any other operation on the 
same machine then will be zero. In such a ease we can 
try for fixing first and second placings of operations on 
this machine and decide the minimum increase in cost to 
be /S^. In this way we can decide -the. minimum increments 
due to third, fourtii end so on placings. If at any stage 
there are r unassigned operations on a machine, ttien for 



deciding the Eininuin increcse in cost due to sll possible 
placing vi'e vdll have to consider rJ pennu tations. We 
would like to fix only the first and the second placings. 
The second placing has been included so that our calcula- 
tions are not sarred if the first placed opera.tion has 
no clashes. It vdll require ttie searching of r(r-l) 
combinations. 


2 

's are calculated for each machine separ- 
ately, and the maximum of them is chosen to be minimum 
increment Z_\ in the total cost. If vm want to consider 
simultaneous effect of on the .late -times then 

we will again have to consiaer r(r-l).ni combinations 
Vvhere n is tiie number of machines. A large number of 
calculations cannot provide complete Scdety against 
backtrack so it is not v.crth it. 


A flow chart for the calculation of Z!^is given 
on the following page. The number of searches for the 
first placing can be cutdovvn to only those operations 
which are in clash vdth an operation having the earliest 
finish time.. 

All operations in the follovving chart refer to the 


unassigned operations. 
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iCi'’-oose e. machine k' 

L : ; — 


Cc'lculpte for all job i, 

Tx = Son of the durations of operations 

on job i that follow the operation, ik 

„ 

Find the operation, iok having | 
the earliest finish time i 
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Then the ^.is chosen to be the maximum of -^2 

for all If machines. The bound will then be cailculated 
as 

m ' m 

1 =11 CJ^ • SJ'i 4 - X ^ (16) 

i=l ~ if=l 

k 4ob-shop problem using the above bounds has 
boen solved in the article 5. of tiiis chapter. 

Before closing this chapter it should be observed 
that th'.e branch and bounc procedure is a powerful tool to 

solve the lob- shop scheduling problem. More than that 
it gives considerable insight into such problems. 



ChiUr'TEfv IV 


i SuBOPTIMAL BRABCh isljp BOOKD PKOCEDOHE 

In the previous chapter ■v«e have seen that to 
eliminate backtraclflng one has to manipulate a la.rge 
amount of data. Even the restricted bounds proposed^ 

the method Vvill reauire considerable computer time for 
large problems. In large problems the effect of schedul- 
ing an operation v.ill give rcdatively small increase in 
in the total cost. The larger the difference between this 
cost and the cost of tiie optimal schedule, the larger the 
b axls'tracking steps ttsat vvould be expecteo. Unfortunately 
even ordinary inciostrial problems are too big for the 
optimal methods proposed in the Isist chapter. This is the 
reason inductries have not shown a keen interest in the 
exact solution of the scheduling problem. Therefore, 
there is a great need of finding some ouick suboptiraal 
procedures. In tiie following aiscussion £■. suboptimal 
procedure is developed on the lines of the branch anu 
bound procedure developed in the previous chapter. 

4.1 S OBOPTIMfL METHOD 

The procedure developed in the previous chapter 
constructs the schedule by proceeding from left to right 
resolving tne clashes accorcing to the lovaest bound 
criteria. Tlois is also tile natursl vvay any scheduler 
will adopt even though he may not be knowing anything 
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about the branch and bound nethod. In place of the 
selection accordin^y to t'ne lowest bound, he often uses 
tbunb rijles or heuristic rule lifee ’first come first 
served’ or ’shortest imminejit operation (SIO)’ rule. The 
heuristically generated schedules are generally one pass 
and no bac^trecMng steps are undertaken. Eovv that we 
have bounding proceuure for any linear cost function we 
can also go for one pass solutions. Therefore the method 
developed in tlse previous chapter, with all backtracM.ng 
steps ignored, can easily replace tlie heuristic and 
thumb rules. 

lilote that in the bounding procedure we have consi- 
dered the effect of next two placings of the unassigned 
operations on all machines. This is the ’look ahead’ 
feature of this procedure. The question that we face 
is, vvhether this is a good guarantee for obtaining a 
near optimad solution. It does not give any such gu0.rantee 
It assures us only of a locally optimum solution. In any 
case its guarantee for a locadly optimum solution is far 
better than by any existing heuristic or thumb rule because 
they do not have ’look ahead’ property. 

4.2 CRITICAL OP ERATIOMS 

Giffler et, el . (l) have found that with the increase 
in tile range of durations of various operations, the 
combinatorial aspects of the problem also increase. This 
is expected because such problems nave higher possibility 
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while resolving clashes among fwo or more 
operations on a machine;, it is quite possible that bounds 
for some nodes have same values. This tie can be resolved 
by selecting an operation which gives larger push to other 
operations in ti'ie schedule. It may also happen that ‘for 
an operation^ whicn gives larger rightward push to other 
operrtions, the bound is only slightly higher than the lowest 
bound. Since tiris operation is lilfely to ease mere clashes 
among other operations it can be a critical operation. The 
selection according to numericailly lowest bound will ignore 
it. A method to take care of such critical operations can 
be suggested as follows ; 

Choose' those nodes Vvhich have bounds within the 
lowest bound ulu.s the i'iicreraent in the 3.ovtfest bound over the 
lovvest bound of the ircocliately preceding branch of tile 
assignment tree. Calculate indices by dividing the in- 
crease of each of 'bhese bounds over -tlrie lowest bound of the 
previous bronch by the displacement that its corresponding 
operation causes to other operations. Ctioose the operation 
corresponding to the lowest index. 

Let us consider -the iob-shop problem as given in the 
Fig.T. Suppose the bounding procedure restricts itself to 
finding the ^ due to placing of only one of the 'next’ 
operations on each machine, as against suggested two placings 
in article 6 of chapter III. The partial tree for -the 
problem will be as follows : 
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The aisplc'cerT'eni: eoeisei hy placing operation 33' in the 
first plrce on the aiochiiie 3 shifts operrtion 13 by 1 
and operation 23 by 5. The totsl di splo cement caused is 
6. Tr e total da splaceesjent caused by operstion 23 is 1, 
Then the inuices a re,. 


ince? for operetion 33 



index fc 


operation z. 


7D - 67 



The index is lov.er for tne operation 33 even though its 
bound value is nur.erically higher than tiiat o-f operation 
23. Note theat o. displacement of 6 by operation 33 does 
not give relatively as higher a increase in tJae bound as 
the operstion 23 givenWith the displaconont of only 
1 unit. Thoreforo we decide on 33 as the critical 
operation and schedule it. Tile resulting schedule after 
till a choice rdves z total cost of 17 which is optiraali 
while that corresponding to the c}:m)ice of operation 23 
givfls a total cost of 20. - The Gnatt charts are as 
follows. 
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Fig. 3 

4. 5 C0MPAKIS0 I4 

A conjperison Ijetvueen the subop timal procedures witic 
and vsitiiout the facility of recogni 2 ,inf the critical operation 
given in Table. 4*1* Th-e earlier procedure is referred to as 
procedure 2 and the later as procedure 1. The bounding 
procedure is the senie as given in article 6 of the./ Chapter III. 

The problems^ selected for comparison are such that 
the differences in the durations are large, which increases 
the combinatorial aspects of the problems. There are two 
types of problems for eacl'i siae^ one of which has identical 
sequences of processing on all Jobs and the other has non- 
identical sequences. The later further increases the com- 
binstorial aspects of the problem. The costs of unit job 
late time and machine idle-time are aliie and unity. The 
problems ere given in Appendix 1. 

The results listed in the table are the costs of 
the overall sqnedules. be can see tiiat the results obtained 
by procedure 1 are in general better than those obtained by 
procedure 2. This shows tiiat the scheme used to detect the 
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criticsl ■opera.tions i-er net su-efi'©S'S^ui » S&a-e- 

is required, in this cLirectiou. 

The tine oy 7044 to solve all tile problems 

in I'fJjle ,4*1 is nearly 480 seconds for bofeh procedures 1 
end 2. Out of tills time nearly 100 seconds are used in 
compiling the program. 

Table 4*1 ciso conperes the results of procedures 
1 and 2 v.itii the one -produced by using heuristic rule SIO 
(giortest Imminent Operation Rule). Accoro.ing to this rule 
all those operations in tiie mix are scheduled which use 
different macliines and are completed in earliest time. 

This rule has beeni found to be very efficient for minimining 
the total time of completion of all jobs(l). Here v/e have 
calculated tr.e costs of botii ooblate'^times end machine 
idle times for the schedule generated by using SIO. 

Vve observe tnat ell tJie results obtained by both 
su'boptimel procedures 1 and 2 are better then that obtained 
by using SIO. In case hhe costs of job latetimes and 
machine idle- times ere nonunity and random, the subop timal 
procedures will be undoubtedly better than the SIO rule. 

Tnis strengtlcens our belief in replacing the heuristic 
rules by the suboptimel procedure 1. 
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Subop time! Suboptimai. 
procedure procedure SIO 

1 2 

1. 

4x3 

identical orders 

53 

56 

56 

2 * 

3x5 

identical orders 

223 

208 

342 

3, 

5x5 

nonicienticel orders 

151 

145 

199 

4 . 

6x6 

identic al 0 r o. e r s 

192 

204 

206 . 

5- 

6 x 6 

nonidentical orders 

215 

220 

247 

6 . 

7x7 

i den ti c ei 0 r cter s 

699 

765 

799 

7* 

7x7 

nonidentical orders 

392 

41 8 

467 

8 . 

8 x 8 

identical orders 

1019 

1007 

1347 

c 

y • 

10 x 10 nonidentical 





orders 

6932 

7 621 

7662 


The cost of schedule 
Table 4-1 


4.5 COMPUTEh, PROGRAM 

A computer program has been Vvritten in FORTfiAM 17 
for tile suboptiraal procedure 1 v,hicii has no facility of 
recognising the critical operations. The program consists 
of one aiain program and three subroutines. The subroutine 
SOB 2 constructs tiie Gnattdiart for any node. The SOB3 
calculates the bounds in the Gnatt chart produced in S0B]_ 
according to the procedure given in article 6 of diapter 
III. The SUB^ finds out an unassigned operation on a laechim 
vvhich has least time of completion with Gnatt chart. This 
is used to find the clashing operaions on a. ma.chine at 
any stage of sclieduling process. 



In the main program the mix is formed and updated 
as the scheduling proceeds. If an operation in the mix 
is not in clash it is directly scheduled v-vithout having 
to call the subroutines. Ihe main program decides upon 
a proper machine in the mix and calls the subroutines to 
construct the Gnatt cherts and calculate the bounds for 
for the clashing operations on the selected machine.! 
decision is taken according to the lowest bound and 
agedn the procedure is repeated. 

Ihe format statements for the input and output of 
data must be provided by the user. Ihe data arrangement 
is as follows : 

(1) Number of jobs and number of machines in 14 format 
on a single card. 

(2) Input format specification for the order matrix. If 
the matrix is 10x6 where' 10 is the number of jobs, 
then a format must be specified as (614) starting 
from the first column cf an IBM card. 

(3) Forma.t specification for the dura.tions matrix. 

If tfce matrix is 10 x 6, then a format like (6F6.1) 
must be specified on a single card starting from 
the first column, Ihe same format is also used for 
the output of the scheduled durations matrix. 

(4) Fbrmat specification for supplying the information 
about the arrival time of jobs and the cost per unit 
jobla tetiraes. Ihe foimf t should be specified as 
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(F6.1, F 0.4} on a single card. , Bie number of 
columns punched should not exceed 24. • The first 
forsict refers to the arrival time of jobs. 

(5) Input data for order matrix. 

(6) Input data for the durations matrix. 

■■(7) Input data for the arrival time of jobs and the 
costs of joblatetines. The order of punching 
these data on a single card must be the same in 
which they are mentioned here. Data for each job 
is punched on separate card. 

(8) Input data for tiie available times on machines and 
the costs of machine idletimes must be punched in 
the same manner as mentioned in (7)., 

The output of the scheduled matrix is arranged in 
the form of a ms.trix, the rows of which are the jobs ana 
the columns indicate tine machines, !]he (i,j)th entry 
in the matrix is the scheduled start time of the operation 
of ith job on jth machine, 

* * 

4.4 ESTlMATlOIj OF T£E OOMPOTER TIME 

In combinatorial problem like the job- shop 
scheduling it is not possible to give a closed form 
relationship betvveen the size of the problem and the 
time required to find the subop timal solution. It all 
depends on how the various operations on different jobs 
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fit v^/ith each other, Ihe uncertainty increases with the 
Increase in range of differences in the processing times, 
we can estimate the time trOken for the worst possible 
case in which we have to resolve clashes among all the 
jobs on all the machines. If there are m jobs and n 
machines, then we have to calculate at worst (m-l)i-l).n 
nodes. 

Bie calculation of nodes in the beginning of the 
schedule requires maximum manipulation of a large data 
while those at the end requires the minimum. The 
average calculations involved in a node which is in the 
middle can be estimated as follows : 

Calculation of SMjj. (k= 1,... n/2) : — DO loops of length 

2 

— nni 

2* ^ — additions. 

4 

Calculation of ^ ^ ^ = l,...n/2) ; ^ *2 

additions. 

Construction of Gnatt chart : ( ~ — 1) shifts or additions. 
The various tests may be expected to take as much time as 
all the above additions will take. Thus a very very crude 
estimation of the equivalent number of additions will be, 

^ ' 4 - : ’ ' ^ 

Generally m is very large compared to n, then this number 
can be approximated as 

, n m^/4 



2 2 


Therefore for all the nod^s t-oproximately numhe:; 

H' 

of equivalent additions are required to be done. 


This number is very strongly dependent upon the 
number of jobs. This number is very crude because it is 
very unlikely that all the unassigned operations 
simultaneously clash over all machines. However, we can 
expect that this maximum estimate should have proportionate 
relationship with the actual time of computation. 


The estimation of minimum time of computation is 
very difficult for general problem in which there is 
nonidentical orders of processing various jobs on 
different machines. For the iaentical-order case 1^-e 
miniraum number will be Cm-l)! nodes. The time required 
for calculrtion cf average node will be nearly same as 
in the v/orst case. But this number is not so informative 
as the maximum number calculated for the worst case, l/hen 
jobs increase from m to m 1, the increase in the cost of 
computation can be found as follows 

ACostoCC. j (m-tl)i - (m)J J 

. , ' ■ ^ , J ■ / . , . ' 

where C is a constant. 

Further simplification gives 


4^ , Cost oC C 





0 2 
n*'m 






that is, Acostos.cost of pre-vious schedule multipli 

by m. 

Ihls gives us the cost of computation as an 
exponential function of the square of the number of 
jobs Vihen the number of machines in the ibop is 
constant. 
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.:NCLu£1( 


’RACT’ 


COATIOi^S 


The optirne-l brench and bound method and the sub- 
cpuii;-?! method, both have linear objective functions. In 
general\the objective function is nonlinear, for example 
the due date problerr. is very common. Likewise tiiere are 
many other constraints like assembly end dis^assembly, 
v'.hich occur almost without fall in eai industrial acti- 
vities. ■ Besides, very often nevv job orders arrive and 
what shruld be the managements reaction twoards it ? All 
these are very practical problems. A few suggestions to 
take account of such situations are given in the follow- 
ing discussion . 

5.1 ARRIViiL OF A giElfe JOB ORDER 

A new job may arrive while 'a batch of jobs is 
still under processing. The management is faced wdth 
the problem vvhether the present batch should be disturbed 
or not# If the new job is a rush order then the .present 
batch have to be disturbed. In this case the old schedule 
has to be scrapped completely and a new branch and bound 
solution must be generated from the arrival time of the 
new job. In case there is planning ahead then one can 
think of disturbing the schedule even prior to the arrival 
■time of the new job. It all depends on the expected profit 
from the new job. From the previous chapter we know that 
the cost of computation is an exponential fun ctioh of the 
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number of Jobs. IJi© iji tha ©o-st o# 

compute t.lon involving e neo job will be. 

L'ost- compiitc.tion = Cost of compu ts,tion of originsl 

schedule multiplied by then 
number of jobs in the previous 
schedule. 

The profit of scheduling a new job can be estimated as 
follows : 

Fit the new job to' the existant schedule without 
disturbing it and calculate tiie cost of the new schedule. 

Let it be Ci. 

Suspend scheduling all the opera.tions from a point 
fro 23 v/hich the new job disturbs the madiine idle- times 
of the previous schedule. Perm the mix of next immediate 
operations at that instant and calculate a node corres- 
ponding to giving priority to tiie new job. Let the 
resulting node value be 

If the cost of the previous schedule be C, then 
C2 - C is the minimum increment in the cost. The difference 
between C2 *■ C and Cj - C is an estimate of the profit from 
from the new job. If this compares favourably with the 
increment in the cost of compu tation, then we should re- 
schedule otherwise not. 



So 


5.2 KiE DATES 

In these ’probieris the objective function is piece- 
wise linear* If the constraints of tiie LPP are linear then 
the region of the solution Vvill still be s convex polyhedron 
or r dimensions where r is the number of variables. Uie 
intersection of the objective function with this polyhedron 
will be a polyhedron of si2,e r or greater dimensions as 
against a polyhedron of r-1 siae in case of linear objective 
function. Integer! nation of any varicible vvill be equivalent, 
to parallel slhifting of tl'.e intersection polyhedron. If 
tile integer takes vilues of 0 or 1 only^then the objective 
function corresponding to 2.ero or one will constitute the 
lower bound for the objective function. Hence the branch 
and bound procedure can still be applied to this case. 

This technique is a very general search technique. 

The bounding procedure given in the third chapter 
can be used for the above problem. The difference will be 
only in selecting the new machine to be considered for 
resolving clashes. A good selection criteria can be the one 
which has the highest derivative of the cost function** 

5.3 ASSMBL2 MD DISASSEMBLY OF JOBS 

The assembly of two or more jobs can be eonsidered 
as yet another job which has its earliest starting time 
dependent on the time of finishing of all the assembled 
jobs. If job i and j are assembled, then for the assembled 
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job .1 ci'.n ^dit ifolxo-v^ijcig -a'anss-t^airit-. ija. ti/e siat. of 

eqii.rtions 



(1) 

2 , 


-t 

(2) 


The constrr:ints vdll be tr-ken care of v/hile constructing 
the Gnatt chart. 

The cost of joblateness is assigned to the assanbled 
job snd not the constituent parts. 

5.4 IBTERROPTIOh OF OPERATIONS 

If Vv’e allow interruption then the number of nodes 
which must be calculated increases. Consider the following 
two jobs in clash over the ktio machine. 

Job 1 
Job 2 

The resolution of clash can be done in the following ways 



The second possibility is a creation of the interrupting 
facility. The bounds are calculated for all these possi» 
bilities and the one with the lowest value of the bound 
is chosen. If the number of jobs in simultaneous clash is 
more than two the possibilities will increase, so also the 
time of computation.. The bounds must also taJce into account 
the interruptibility of operations while calculating SE&jf's 
■and the incremntal costs. 




5.5 MORE ThAIM OKE OF A TCINi^ 

Biis is a very practical situation to occur in any 
type of industry.There may be more than one lathe of the 
sv.me kind in a. workshop. Just like the case of interrupted 
operations it also increases the number of nodes that must 
be calculated in order to find the lowest bound, 

liae node calculations will be slightly different 
than that given in Chapter III. The SMjj*s will now be 
calculated as the minimum time of finishing the last 
operation on the group of machines of kind.k, taking into 
account all the possibilities of simultaneous operation 
of various jobs on these machines. Since this will require 
consider search time it will be advisable to consider the 
time of finishing the last operation on the kth machine in 
the present partial schedule. 

The cfilculation of incremental costs due to placing 
of next two unassigned jobs on the kth machine will also 
become cumbersome. The incremental costs thus obtained 
will also be of veiylittle help because a lot of compact- 
ness in the schedule is to be expected. Therefore, the 
cBlculation of bound becomes very simple. 

5.6 ALTERNATE ROUTES 

Some jobs may have alternate routes of processing 
on various machines. Note that in the fixed-route case 
the SJi’s, the time of finishing the last operation on the 
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ith- 3 ob, if; fixed. In tiiis case the KLlX T^e 

calculated to be the minimum time of finishing among 
the alternate routings of the unassigned operations on the 
ith ;]ob. The various combinations of routes will complicate 
the calculation of SSdjf and also the incremental costs due to 
inter- job interactions. Therefore, as suggested in the 
article 5.5, the calculation of incremental costs is suspended 
altogether and the are chosen to be the time of finishing 
the last scheduled operation on the kth machine plus the sum 
of the processing times cf the remaining operations on this 
machine. 

5.7 C0I3CLUS1OM 

The search procedure given in the Chapters III 
and IV can incorporate many practical situations with 
very little change. 



APPEKL-IX I 
SAMPLE PROBLMS 


1. 

1(3) 

2(10) 

3(4) 





1(4) 

2(1} 

3(5) 





1 (4) 

2(6) 

•3(6) 





1(6) 

2(4) 

3(1) 




2. ' 

1(3) 

2(6) 

3(9) 

4(1) 

5(8) 



1(5) 

2(16) 

3(6) 

4(^) 

5(10) 



1(9) 

2(2) 

3(7) 

4(10) 

5(15) 



1(3) 

2(9)- , 

3(14) 

4(2) 

5(15) 



1(1) 

2(8) 

3(25) 

4(22), 

5(3) 



3(3) 

2(6) . 

, 5(9) 

1(1) 

4(8) 



4(5) 

1(16) 

2(6) 

3(8) 

5(10) 



5(9) 

2(2) 

4(7) 

1(10) 

3.(15) 



4(3) 

1(9) 

2(14) 

3(2) 

5(15) 



2(1) 

3(8) 

4(25) 

5(22) 

1(3) 


4. 

1(1) 

2(3) 

3(6) 

4(7) 

5(3) 

6(6) 


1(8) 

2 ( 5 ) 

3(10) 

4(10) 

5(10) 

6(4) 


1(5) 

2(4) 

3(8) 

4(9) 

5(1) 

6(7) 


1(5) 

2(5) 

3(5) 

4(3) 

5(8) 

6(9) 


1(9) 

2^3) 

3(5) 

4 '(4) 

5(3) 

6(1) 


1(3) 

2(3) 

3(9) 

4(10) 

5(4) 

6(1) 



1-2 


3(1) 

1(3) 

2(6) 

4(7) 

6(3) 

5(6) 


2(8) 

3(5) 

5(10) 

6(10) 

1 (10) 

4(4) 


3(5) 

4(4) 

6(8) 

1(9) 

20.) 

5(7) 


2(5) 

1(5) 

3(5) 

4(3) 

5(8) 

6(9) 


3(9) 

2(3) 

5(5) 

6(4) 

1(3) 

4(1) 


2(3) 

4(3) 

6(9) 

i(ic) 

5(4) 

3(1) 


1(12) 

2(3) 

3(15) 

4(9) 

5V5) 

6(6) 

7(2) 

1 (3) 

2(6) 

3(21) 

4(8) 

5(4) 

6(12) 

7(10) 

1(8) 

2(28) 

3(6) 

4(40) 

5(2) 

6(8) 

7(9) 

1(7) 

2(25) 

3(14) 

4(6) 

5(9-> 

6(1) 

7(13) 

1 (26) 

2(9) 

3(3) 

4(12) 

5(5) 

6V24) 

7(7) 

1(5) 

2(12) 

3(6) 

4U6) 

5(9) 

6(34) 

7(7) 

1(12) 

2(16) 

3(8) 

4(7) 

5(19) 

6(6) 

7(3) 

1 (12) 

2(3) 

3(15) 

4(9) 

5(5) 

6(6) 

7(2) 

3(3) 

2(6) 

6(21) 

5(8) 

1(4) 

4(12) 

7(10) 

4(8) 

1(28) 

7(6) 

6(40) 

2(2) 

3(8) 

5(9) 

5(7) 

2(25) 

3(14) 

4(6) 

6(9) 

7(1) 

1(13) 

4(26) 

1(9) 

6(3) 

7(12) 

3(5) 

2(24) 

5(7) 

2(5) 

7(12) 

4(6) 

5(16) 

1(9) 

3'^34) 

6(7) 


7(12) 6(16) 5(8) 4(7) 2(19) 3(6) 1(3) 




g. 1(2) 2(13)' 3(24) 4(10) 5(5) '6(6) 7(9) S(10) 

1(36) 2(6) 3(15) 4(5) 5(9) 6(3) 7(8) 8(24) 

1(5) 2(25) 3(23) 4(6) 5(38) 6(9) 7(3) 8(2) 

1(9) 2(12) 3^3) 4(35) 5(6) 6(9) 7(1) 8(3) 

1(5) 2':'^ 3(19) 4(1) 5(5) 6(12) 7(15) 8(3) 

1(9) 2(13) 3(5) 4(18) 5(9) 6(17) 7(20) 8(30) 

1(5) 2(3) 3(2) 4(25) 5(36) 6(9) 7(1) 8(7) 

1(6) 2(9) 5(19) 4(4) 5(22) 6(6) 7(7) 8(28) 


9. 1(29) 2(78} 3(9) 4(36) 5(49) 6(11) 7(62) 8(56) 9'(44) 10(21) 

1(43) 3(90) 5(75) lO(IL) 4(69) 2(28 ) 7(46) 6(46) 8(72) 9(30) 

2(91) 1(85) 4(39) 3(74) 9(9C) 6(10) 8(12) 7(89) 10(45) 5(33) 

2(81) 3(95) 1(71) 5(99) 7(9) 9(52) 8(85) 4(98)10(22) 6(43) 
3(14) 1(6) 2(22) 6(61) 4(26) 5(69) 9(21) 8(49)10(72) 7(53) 
3(84) 1 .( 2 ) 6(52) 4(95) 9(48)10(72) 1(47) 7(65) 5(6) 8(25) 

2(46) 1(37) 4(61) 3(13) 7(32) 6(21)10(32) 9(99) 8(30) 5(55) 

3(31) 1 ( 86 ) -2(40) 6(74) 5(32) 7(88) 9(19) 10(48) 8(36) 4(79) 

1(76) 2(69) 4(76) 6(51) 3(65)00(11) 7(40 ) 8(89) 5(26) 9(74) 

2(85) 1(13) 3(61) 7(7) 9(64)10(76) 6(47) 4'^52) 5(90) 8(45) 
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